Neutrino mass limits from SDSS, 2dFGRS and WMAP 
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We investigate whetlier cosmological data suggest the need for massive neutrinos. We employ 
galaxy power spectrum measurements from the Sloan Digital Sky Survey (SDSS) and the Two 
Degree Field Galaxy Redshift Survey (2dFGRS), along with cosmic microwave background (CMB) 
data from the Wilkinson Microwave Anisotropy Probe (WMAP) and 27 other CMB experiments. 
We also use the measurement of the Hubble parameter from the Hubble Space Telescope (HST) 
Key Project. We find the sum of the neutrino masses to be smaller than 0.75 eV at 2a (1.1 eV at 
3a). 
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Neutrino oscillation experiments provide substantial 
evidence that the three known neutrinos have a com- 
bined mass (S = ^ m^) of at least about y^Jmf ~ 0.045 
eV^, but are completely insensitive to the eigenvalue of 
the lightest neutrino mass eigenstate [1]. The neutrino 
mass spectrum may be hierarchical (mi ^ m2 <C m^, 
TO3 <^ mi < TO2) or quasi-degenerate (mi ~ m2 — ma), 
depending on whether the lightest eigenmass is close to 

or 3> \/ Srn^, respectively. The nature of the spectrum 
is important to neutrino mass model-building, the con- 
tribution of neutrinos to dark matter, and the viability 
of observing neutrinoless double beta-decay if neutrinos 
are Majorana [2]. 

Major progress towards the determination of the abso- 
lute neutrino mass scale has been made with laboratory 
based measurements and with post- WMAP cosmologi- 
cal data. Tritium beta decay experiments constrain E 
to be smaller than 6.6 eV at 2a [3] which can be im- 
proved to S « 1 eV by the KATRIN experiment [4]. 
Large scale structure (LSS) data from SDSS [5] combined 
with CMB data from WMAP [6] alone, yield E < 1.7 
eV at the 95% C. L. [7], with no strong priors or as- 
sumptions. With 2dFGRS data [8], CMB data and sig- 
nificantly stronger priors, a 95% C. L. upper limit of 

1 eV was found in Ref. [9]. An even stronger bound, 
S < 0.63 eV was obtained by the WMAP collabora- 
tion from a combination of 2dFGRS and CMB data and 
the 2dFGRS measurement of the galaxy bias parameter 
b = Pg{k) / Pjn{k) , where Pg and Pm are the galaxy and 
matter power spectra, respectively. The 2dFGRS col- 
laboration measured a scale-independent bias over scales 
A: ~ 0.1 - 0.5 /i/Mpc [10]; the WMAP collaboration 
adopted this bias for k < 0.2 /i/Mpc. A preference for 
massive neutrinos found in Ref. [11] is controversial be- 
cause of the input of a linear clustering amplitude erg 
(defined as the rms mass fluctuations in spheres of ra- 
dius 8 h^^ Mpc) [12]; at present there is no experimental 



consensus on the determination of erg {e.g., see Table 5 
of Ref. [7] ) . All the above cosmological constraints were 
placed under the assumption of a flat Universe in accord 
with the predictions of inflation. 

We analyze a large set of LSS and CMB data at 
scales where the matter power spectrum is linear, i.e., 
k < 0.15 /i/Mpc. We include the power spectrum de- 
terminations from 205,443 and 147,024 galaxy redshifts 
measured by SDSS and 2dFGRS, respectively. The CMB 
data comprise all WMAP data and a combination of 
the 151 band power measurements from 27 other CMB 
experiments [13] including CBI [14] and ACBAR [15], 
with multipoles / up to 1700 (or k ~ 0.15 h/Mpc [16]). 
Throughout, we impose a top- hat prior on the Hubble 
constant h (Ha = lOOh km/s/Mpc), from the HST [17]. 
We do not include Ly-a forest data [18] in our analy- 
sis because an inversion from the flux power spectrum 
to the linear power spectrum is nonlinear and model- 
dependent [19]. 

EfTects of neutrino mass on the power spectrum: 

Neutrinos of eV masses are relativistic when they de- 
couple, and so their final number density is independent 
of their mass, rii, = 3/lln^. Since (E^) = 2.7Tj, and 
il.yh'^ is essentially the energy density of the CMB with 
T-y — 2.725 K, Ury is known, and 



94.1 eV 



(1) 



^ Sm^ is the mass-squared difference of atmospheric neutrino oscil- 
lations. The mass-squared difference of solar neutrino oscillations 
is significantly smaller with yj (5m| ~ 0.008 eV. 



Neutrinos freestream on scales smaller than their Jeans 
length scale, which is known as the freestreaming scale. 
While neutrinos freestream, their density perturbations 
are damped, and simultaneously the perturbations of 
cold dark matter and baryons grow more slowly because 
of the missing gravitational contribution from neutrinos. 
The freestreaming scale of relativistic neutrinos grows 
with the horizon. When the neutrinos become nonrel- 
ativistic, their freestreaming scale shrinks, they fall back 
into the potential wells, and the neutrino density per- 
turbation resumes to trace those of the other species. 
Freestreaming suppresses the power spectrum on scales 
smaller than the horizon when the neutrinos become non- 
relativistic. (For eV neutrinos, this is the horizon at 



2 



matter-radiation equality). Lighter neutrinos freestream 
out of larger scales and cause the power spectrum sup- 
pression to begin at smaller wavenumbers [20] , 

fc„.=. 0.026 (^^y\pc-i, (2) 

assuming almost degenerate neutrinos. Here, lom = 
Q.Mh'^ is the total matter density (which is comprised 
of baryons, cold dark matter and massive neutrinos). On 
the other hand, heavier neutrinos constitute a larger frac- 
tion of the matter budget and suppress power on smaller 
scales more strongly than lighter neutrinos [21]: 

where = 17^/17 a/ is the fractional contribution of neu- 
trinos to the total matter density. 

Analyses of CMB data are not sensitive to neutrino 
masses due to the fact that at the epoch of last scatter- 
ing, eV mass neutrinos behave essentially like cold dark 
matter. (WMAP data alone allow the dark matter to be 
entirely constituted by massive neutrinos [7]). However, 
an important role of CMB data is to constrain other pa- 
rameters that are degenerate with S. Also, since there 
is a range of scales common to the CMB and LSS ex- 
periments, CMB data provides an important constraint 
on the bias parameters. Sensitivity to neutrino masses 
results from the complementarity of galaxy surveys and 
CMB experiments. 

Figure 1 shows that the suppression of power caused by 
massive neutrinos is much greater for the galaxy power 
spectrum than for the CMB TT spectrum. We do not 
show the effect of neutrino masses on the CMB TE spec- 
trum because it is tiny. Note that we have normalized 
the spectra to emphasize the power suppression at small 
scales. 

Analysis: 

We compute the CMB TT and TE power spectra 
5Ti = 1(1 + l)Ci/27r, and the matter power spectrum 
Pm{k), all in the linear approximation, using the Code for 
Anisotropics in the Microwave Background or CAMB [22] 
(which is a parallelized version of CMBFAST [23]). We 
assume the Universe to be flat, fttot = 1, that the dark 
energy is in the form of a cosmological constant A, and 
that there are three neutrino species. We calculate the 
angular power spectra on a grid defined by ujm, fv, the 
baryon density lob = flsh^, the Hubble constant h, the 
reionization optical depth t, and the spectral index Ug of 
the primordial power spectrum. 

We employ the following grid: 

• 0.05 < ujm < 0.27 in steps of size 0.02, and ujm = 
0.14, 0.16. 

• < /i/ < 0.15 in steps of size 0.01. 
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FIG. 1: Upper panel: CMB TT power spectra. Lower 
panel: Galaxy power spectra. The curves are the spec- 
tra for E — 0.28 eV (solid, best-fit parameters; the galaxy 
power spectrum is shown for the SDSS best-fit normaliza- 
tion), E = 1.5 eV (dotted) and S = 3 eV (dashed). The latter 
two spectra have all other parameters (except the normaliza- 
tion, bias parameters and flcdm) fixed at the best-fit values. 
All curves are for a flat universe. The CMB TT spectra are 
normalized to have identical powers at the first peak. The 
galaxy power spectra are normalized to have identical powers 
at fc = 0.017 h/Mpc. In the upper panel, the data points 
marked by circles (squares) represent the binned TT spec- 
trum from WMAP (pre- WMAP experiments). In the lower 
panel, the data points marked by circles (squares) represent 
the galaxy power spectra from the 17 SDSS (32 2dFGRS) 
bands used in our analysis. 



• 0.018 <UJB < 0.028 in steps of size 0.001. 

• 0.64 <h< 0.80 in steps of size 0.02. 

• < T < 0.3 in steps of size 0.025. 

• 0.8 < Us < 1.2 in steps of size 0.02. 
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• The normalization of the primordial power spec- 
trum Ag, is a continuous parameter. 

• The bias parameters bsoss and b2dF are scale- 
independent^ and continuous. 

The siippression of small scale power depends directly 
on fi, and indirectly on S. From Eq. (3), ojm is strongly 
degenerate with S, requiring independent knowledge of 
UJM to break the degeneracy. The SDSS collaboration 
only used WMAP data to provide this information and 
found the somewhat weak, but conservative, 95% C. L. 
bound S < 1.7 eV [7]. Their analysis also gave a la- 
constraint h = 0.645j;o;o40) which lies at the lower end 
of the HST measurement h = 0.72 ± 0.08 [17]. It is 
known that less stringent constraints on E are obtained 
for lower values of h [25] because CMB data then allow 
larger Hm [26]^. Aside from the fact that we are using a 
significantly larger dataset than the SDSS collaboration, 
we expect our analysis to yield a stronger constraint on S 
simply because we constrain h by the HST measurement. 
Note that our grid-range for lum is almost identical to a 
combination of the 3a range of allowed by SN la 
redshift data [28], and the HST prior on h. 

In our analysis, we conservatively include only the 
first 17 SDSS band powers, for which 0.016 < k < 
0.154 /i/Mpc, and the power spectrum is in the linear 
regime. We use the window functions and likelihood code 
provided by the SDSS collaboration [5] , and leave the bias 
parameter bsDSS free. 

For 2dFGRS, only 32 band powers with 0.022 < < 
0.147 h/Mpc are included. The window functions and 
covariance matrix have been made publicly available by 
the 2dFGRS collaboration [8]. The bias parameter b2dF 
is left free. 

The WMAP data are in the form of 899 measurements 
of the TT power spectrum from ? = 2 to ? = 900 [29] 
and 449 data points of the TE power spectrum [30] . We 
compute the likelihood of each model of our grid using 
Version 1.1 of the code provided by the collaboration [31]. 
The WMAP code computes the full covariance matrix 
under the assumption that the off-diagonal terms are sub- 
dominant. This approximation breaks down for unreal- 
istically small amplitudes. When the height of the first 
peak is below 5000 fxK^ (which is many standard devia- 



^ On large scales, galaxy bias is expected to be a scale-independent 
constant [24]. This has been confirmed by the 2dFGRS collabo- 
ration [10]. 

^ As a proof of principle, it was shown in Rcf. [25] that a flat non- 
ACDM model with h = 0.45, = 1 and S = 3.8 eV, provides 
as good a representation of the 2dFGRS and pre-WMAP CMB 
data as a flat ACDM model with h = 0.7, Qm = 0.3 and masslcss 
neutrinos. This is true even including the WMAP data provided 
different power-laws are used to describe the spectrum for I above 
and below the first peak [27]. 



tions away from the data) , only the diagonal terms of the 
covariance matrix are used to compute the likelihood. 

We include the combined CMB data from pre-WMAP 
experiments, by using the 28 pre-WMAP band powers, 
the window functions and the correlation matrix com- 
piled in Ref. [13]. 

We obtain the x-a range of E, by selecting those pa- 
rameter sets with Ax^ = ~ Xmin — ; after minimiz- 
ing over all other parameters. 

Results and conclusions: 
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FIG. 2: Ax^ vs E from an analysis of SDSS, 2dFRGS, WMAP 
and other CMB data with bsoss and b2dF free. 

Our analysis yields the best-fit parameters fi, = 0.02, 
LVM = 0.15, LUB = 0.023, h = 0.66, r = 0.075, = 0.96, 
Ag = 21.38 X 10-1° with = 1499.83 for 1425 - 9 = 
1416 degrees of freedom. From the normalizations of the 
power spectra required to fit the SDSS and 2dFGRS data, 
we obtain bsoss = l-13to.i3 and b2dF = l-ZOto.ls (lo" 
ranges) . 

Figure 2 shows Ax^ versus E. The neutrino mass 
bound is E < 0.75 eV at 2a (E < 1.1 eV at 3a). Thus, 
cosmological data do not require a significant neutrino 
dark matter component, and are increasingly rejecting a 
quasi-degenerate neutrino mass spectrum. 

Eventually, lensing measurements of galaxies and the 
CMB by large scale structure are expected to probe a hi- 
erarchical neutrino mass spectrum with S w 0.04 eV [32]. 
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